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The structural phase behavior of polymer brushes, single-component linear homopolymers grafted onto a pla- 
nar substrate, is studied using the molecular Monte Carlo method in 3 dimensions. When simulation parameters 
of the system are set in regions of macrophase separation of solution for the corresponding non-grafted ho- 
mopolymers, the grafted polymers also prefer segregation. However, macrophase separation is disallowed due 
to the spatially-fixed grafting points of the polymers. Such constraints on the grafting are similar to connecting 
points between blocks of non-grafted diblock copolymers at the microphase separation in the melt state. This 
results in "microphase separation" of the homopolymer brush in the lateral direction of the substrate. Here we 
extensively search the parameter space and reveal various lateral domain patterns that are similar to those found 
in diblock copolymer melts at microphase separation. 



Polymer chains grafted onto substrates, i.e. polymer 
brushes, have been extensively studied in physics, chemistry, 
materials science, and other scientific and industrial fields. 
Polymer brushes are widely utilized for applications such as 
wetting, adhesion, surface patterning, and colloidal stabiliza- 
tion. Recently, sophisticated polymer brushes, such as brushes 
composed of diblock copolymers and of binary mixture of ho- 
mopolymers, have drawn significant attention [1, 2]. Poly- 
mer brushes on spherical substrates, which are referred to as 
polymer-grafted colloidal particles, have also attracted wide 
interest in various scientific and industrial fields [3-5]. 

The phase behavior of these complicated brushes is sig- 
nificantly dependent on the molecular architecture, shape of 
the substrate, ratios between the numbers of polymers of each 
polymer species, and other characteristics of each system. In 
the present work, the most basic and simplest polymer brush, 
i.e. single-component linear homopolymers homogeneously 
grafted onto a planar substrate, is simulated using the molecu- 
lar Monte Carlo method to determine the universal and struc- 
tural phase behavior of polymer brushes. 

The phase behavior of a solution of single-component non- 
grafted homopolymers underlies the phase behavior of the 
present brush. When the chemical and physical parameters 
of this polymer solution are set in regions within binodal lines 
of the phase diagram, the non-grafted polymers are segregated 
into a low-density domain and a high-density domain, which 
indicates macrophase separation. When the parameters of a 
layer of the grafted polymers are set in these regions within the 
binodal lines of the non-grafted polymers, the grafted poly- 
mers also prefer segregation. However, due to a topological 
constraint of the spatially-fixed grafting points of each poly- 
mer, macrophase separation of the grafted polymers is disal- 
lowed. This constraint has similarity to the topological con- 
straint of connecting points between blocks of non-grafted di- 
block copolymers at phase separation in the melt state. These 



non-grafted diblock copolymers exhibit microphase separa- 
tion and various domain patterns [6]; therefore, the grafted 
single-component homopolymers also exhibit phase separa- 
tion in the lateral direction of the substrate, which results in 
a variety of 2-dimensional domain patterns of the grafted ho- 
mopolymers similar to the patterns observed with microphase 
separation of the non-grafted diblock copolymer melts. These 
results indicate that the binodal lines of the grafted homopoly- 
mer system are changed from those of the non-grafted ho- 
mopolymers. In the present work, we extensively search the 
parameter space and reveal the various domain patterns of 
a 2-dimensional "microphase separation" of the grafted ho- 
mopolymers by measuring the lateral density distribution of 
the brush. 

Early research on single-component linear homopolymers 
grafted onto planar substrates predicted and confirmed one 
domain pattern; small circular clusters of the grafted poly- 
mers [7-13]. These clusters are distributed on the substrate 
and are referred to as pinned micelles. The lateral phase sep- 
aration of the brush is termed microphase separation. How- 
ever, determination of the spatial arrangement of these pinned 
micelles was left for future study. Moreover, in this early re- 
search, the possibility of the other domain patterns was missed 
and absent. Here, we reveal a variety of domain patterns and 
complete the structural phase diagram by extensive simula- 
tions. 

In the present work, we employ a 3-dimensional solvent- 
free coarse-grained model [14, 15] of single-component linear 
homopolymers proposed by Miiller and Daoulas [16, 17]. In 
this model, explicit solvent molecules are integrated out and 
replaced with an effective non-bonded interaction between 
solute molecules. This drastically reduces the degrees of 
freedom of the system and significantly diminishes the com- 
putational requirements for simulation of the polymer brush 
placed in 3 dimensions. The linear molecular architecture of 
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the homopolymers is described by a bead-spring potential: 
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where £ B r denotes the thermal energy, r,(s) represents the co- 
ordinate of the i-th coarse-grained segment of the 2-th poly- 
mer, and N is the number of segments per polymer. R e , which 
denotes the root mean square of the end to end distance of an 
ideal chain with the same molecular architecture, is chosen as 
the unit length, //bond represents the bonded interaction poten- 
tial between the segments. On the other hand, the free energy 
of non-bonded interactions in the solvent-free model is given 
as a functional of the local segment density. A third-order ex- 
pansion of the non-bonded interaction free energy in a form 
of powers of the local segment density is employed: 
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where p s (r) is the local volumetric number density of the 
segments at the spatial position r. The positive constants, v 
and w, correspond to the attractive and repulsive interaction 
strengths among the segments, respectively. In the follow- 
ing calculation, dimensionless physical quantities are utilized: 
the local polymer density p' p (r) = p s (r)Rl/N, the parameter 
W - wN 3 /R e 6 , and the parameter v' = vN 2 IR? e . These dimen- 
sionless quantities reduce eq. (2) to: 
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First, we discuss the phase behavior of a solution of the 
non-grafted polymers modeled using the solvent-free system. 
The dimensionless average volumetric polymer density is de- 
noted by p' p = tipRllV, where n p and V are the number of 
polymers and the system volume, respectively. At extremely 
large v' and finite w', the polymers aggregate, which results 
in one dense droplet floating in a dilute gas, which indicates 
macrophase separation of the polymers. In contrast, at finite 
v' and extremely high w', the system is a homogeneous phase. 
The binodal lines in the p' V -plane at fixed W are calculated 
based on mean-field theory [16, 18]. The critical point of this 
phase diagram is denoted by p' p - p' pc and v' = v' c [16]: 
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At high average polymer density, the fluctuation of the poly- 
mer density becomes small, so that the mean-field theory pro- 
vides an accurate prediction of the binodal lines. As an ex- 
ample, the result of the mean-field calculation at fixed W = 
0.0001 is shown in Fig. 1. At this value of w', p' pc = 70.71. 
This high p' pc [16] provides a good prediction of the bin- 
odal line. Based on this phase diagram, the polymer brush 
is simulated and the structural phase diagram is constructed at 
W = 0.0001. 
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FIG. 1. Binodal lines of a system composed of non-grafted poly- 
mers at fixed w' = 0.0001, constructed based on the mean-field the- 
ory. The critical point, (p' pc , v' c ), is denoted by a circle. The broken 
line represents the asymptotic behavior of the binodal line at high 
polymer density. The validity of this phase diagram is confirmed by 
molecular Monte Carlo simulation. 



In this simulation, p' p (r) is calculated via a collocation lat- 
tice. For this purpose, the system box is partitioned into a 
cubic lattice, denoted by {a}, with AL of the grid spacing. 
According to related particle-to-mesh methods in electrostat- 
ics [19-21], the densities at each grid point, a, are determined 
as: 
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The function in this equation, cD (r, a), allocates the particles 
to grid points, according to a linear assignment function, 



O (r, a) = ]~[ g (d a ) , 



(6) 



1 - \d a \/AL for \d a \ < AL 
g{d a ) = \ . (7) 

otherwise 



where d a - r a -a a denotes the distance between the grid point, 
a, and the segment position, r, along the Cartesian direction, 
a. The grid spacing, AL, defines the range of the non-bonded 
interaction. Utilizing this grid-based density, the non-bonded 
interaction of eq. (3) from p' p {a) is calculated by substituting 
the summation over the lattice nodes Y^d^L 3 for the integral 
J y dV. This calculation of the non-bonded interaction via the 
collocation grid is computationally advantageous, because a 
segment in this simulation interacts with many neighbors. 

Monte Carlo simulations are performed with the canonical 
ensemble in 3 dimensions via the standard Metropolis algo- 
rithm [22, 23] at the fixed w' = 0.0001. This value of W 
is the same as that for the phase diagram in Fig. 1. kgT 
is chosen as the unit energy. (L x ,L y ,L z ) denotes the size 
of the rectangular parallelepiped system box and is fixed at 
(L x = l2.0R e ,Ly = l2.0R e ,L z = 6.0R e ). Aperiodic boundary 
condition is applied to the system. The system box is placed 
in spatial regions of < a < L a . N — 32 and AL = (l/6)R e 
are fixed. The Mersenne Twister algorithm is selected as a 
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random number generator for the simulations [24-26]. In 
one simulation step, a particle is picked at random and given 
a uniform random trial displacement within a cube of edge 
length 2AL. A Monte Carlo step (MCS) is defined as n p N 
trial moves, during which each particle is selected for the trial 
displacement once on average. 

One end segment of each polymer is grafted onto and fixed 
at a homogeneous square lattice on a thin hard planar wall 
placed at z = 0. This means that the polymers are homo- 
geneously grafted onto the xy-plane of the system box. This 
substrate disallows the polymers to pass through. No other 
interactions, such as coulombic and frictional interactions, act 
between the substrate and the polymers. The other segments 
of the grafted polymers are built in z > and form a poly- 
mer layer in an interval of < z < lR e - The average vol- 
umetric polymer density in this region of the layer, defined 
as p p = n p /(L x L y x lR e ), corresponds to the average volu- 
metric polymer density in the system of the non-grafted poly- 
mers. Simulations are performed for various value sets of the 
parameters (p' p - p p Rl,v') of the brush, and the structural 
phase diagram is constructed in the p' p v' -plane. We have con- 
firmed that this phase diagram for the brush is not significantly 
changed when the graft points of the polymers are randomly 
distributed over the substrate using uniform random numbers, 
which results in an inhomogeneous spatial distribution of the 
graft points. 

In the initial state, the conformation of each polymer is ran- 
domly arranged according to a Gaussian distribution, which 
results in random coils [27] with &lR e of the root mean square 
end-to-end distance. 

To characterize the domain patterns, the local lateral density 
of the grafted polymers is defined as: 
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which has a dimension of areal density. The average areal 
density of the polymer brush is denoted by pjf' = n p /L x L y . 
Due to the constraint of the graft, except at extremely high 
V , the periodicity of the domain patterns is ~R e , because the 
bonded interaction energy and conformational entropy of the 
polymers disallow the polymers to stretch ?s>R e or to shrink 
^R e . 

Here, we discuss the simulation results for the brush. The 
phase diagram of the brush, which is traced on the phase dia- 
gram of the non-grafted polymers given in Fig. 1, is shown in 
Fig. 2. In regions of v' /v' c smaller than a chain line drawn 
in this phase diagram of the p^,v'-plane, the homogeneous 
phase of the brush is maintained until the simulation halts at 
~1 x 10 6 MCS. As an example of this homogeneous structure, 
P P \x,y)/p p ' at p' p lp' pc = 1-02 and v' ' [V c = 0.5 is presented in 
Fig. 3. On the other hand, in regions of v' /v' c larger than the 
chain line, p„ (x,y) indicates the 2-dimensional microphase 
separation of the grafted polymers. The domains are sponta- 
neously created at less than ~4x 10 5 MCS and remain until the 
simulation halts at ~1 x 10 6 MCS, without significant change 
of the positions, shapes, and the sizes. This indicates that the 




FIG. 2. Phase diagram of the polymer brush, traced on the phase 
diagram of the non-grafted polymers given in Fig. 1. A chain line 
denotes the binodal line of the brush. White squares represent the 
homogeneous phase of the brush. Diamonds, triangles, and crosses 
represent the hexagonal, lamella, and inverse hexagonal structures, 
respectively. These three structures are simultaneously observed in 
the system at the point of the grey square. The grey circle denotes 
the inverse nematic structure and the asterisk is a structure possibly 
similar to, but different from, the lamellae. 



12 



Q; 



10 - 
8 - 
6 - 
4 - 
2 - 
- 



I 1.8 

- 1.6 

- 1.4 
■ 1.2 

- 1 

- 0.8 
I 0.6 







4 6 

x/R c 



10 12 



FIG. 3. Lateral polymer density distribution, pp i \x,y)/p P i \ at 
p' p /p' pc = 1.02, VIV C = 0.5, and 1 x 10 6 MCS. n p = 10404. 



chain line denotes the binodal line of the polymer brush. This 
binodal line is shifted from that of the non-grafted polymers 
toward large v' jv' c . However, these two binodal lines are still 
qualitatively consistent, as expected. Next, we discuss the do- 
main patterns for the 2-dimensional microphase separation of 
the polymers. 

At p' p lp' pc = 1.02, circular clusters assembled in the hexag- 
onal structure are observed at v'/v' c = 2.0, as shown in 
Fig. 4(a). With this value of p' p /p' pc , the number of poly- 
mers becomes n p = 10404 in the present simulation sys- 
tem. This hexagonal structure is also observed in the region of 
0.509 < p' p lp' pc < 1.98 and 2.0 < v'/v' c < 5.0, which is repre- 
sented by diamonds in Fig. 2. At high v' /v' c in this region, rear- 
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FIG. 4. p { p\x,y)/p i p ) at 1 x 10 6 MCS. Results at (a): (p' p /p' pc = 1.02, v'/v' c = 2.0), (b): (1.02,3.0), (c): (1.98,3.0), (d): (3.92,5.0), (e): 
(2.97, 3.0), and (f): (3.62,4.0). 



rangement of the clusters requires a long time, due to the high 
interfacial energy barrier. This results in long-lived defects in 
the hexagonal structure and inhomogeneity of the cluster size, 
as shown in Fig. 4(b), as an example, at p' p /p' pc = 1-02 and 
v'/v; = 3.0. 

When p p lp'p C is raised from the region of the hexagonal 
structure, the clusters increase in size and merge with each 
other. However, due to the constraint of the graft, the period- 
icity of the domain patterns is constrained to ~R e . Therefore, 
large circular clusters with sizes arranged in the hexag- 
onal structure are not created in the system. Therefore, when 
P'p/Ppc rises, the hexagonal structure changes to another struc- 
ture, which is denoted by triangles in Fig. 2. Figure 4(c) gives 
an example of this structure, i.e. 2-dimensional lamella struc- 
ture, which is sampled at p' p /p' pc - 1.98 and V /v' c = 3.0. 

At V /v' c = 5.0, the inverse hexagonal structure is deter- 
mined in the interval of 2.97 < p' p /p' pc < 3.92. In this struc- 
ture, circular voids are created in the dense polymer layer of 
the brush, as shown in Fig. 4(d). However, at p' p /p' pc = 2.51, 
the hexagonal, lamella, and inverse hexagonal structures si- 
multaneously appear in the system. 

At p' p lp' pc = 2.97 and v'/v' c = 3.0, sausage-like voids are 
formed in the dense polymer layer, as shown in Fig. 4(e). 
These anisotropic shaped voids lie in the same direction lo- 
cally, which is similar to the nematic phase of liquid crystals. 

At p' p lp' pc - 3.62 and v' /v' c = 4.0, a structure similar to 
the lamellae is observed. However, in this structure, low- 



density layers of the lamellae compose junction points similar 
to tripods, as presented in Fig. 4(f). Thus, each junction point 
is composed of and homogeneously connected with the 3 low- 
density layers, which indicates that a structure different from 
the lamellae is evident at this state point. 

Finally, the angular distribution of the end-to-end vectors of 
the grafted polymers is examined. When the graft point of a 
polymer is located in regions of low p p (x, y) in 2-dimensional 
microphase separation, this polymer is laterally stretched to- 
ward regions of high p^ix,^). To confirm this effect, the lati- 
tude 9, of the end-to-end vector, i.e. a vector from the grafted 
segment to the free end of each polymer, is measured. In this 
simulation, the latitudinal distribution defined as, 

,™ ._ Platitude (0) d0 

1 {cos 9 - cos (6 + d9)} ' 

is calculated, where platitude (#) denotes the probability density 
of 9. When the vectors are homogeneously distributed over all 
directions for the latitude and longitude, f(ff) becomes homo- 
geneous, i.e. f(ff) = 1 over the entire interval of < 9 < n. 
Therefore, the inhomogeneity of the distribution in the latitu- 
dinal direction is quickly evident in f(ff). 

Figure 5 shows f(6) at p' p /p' pc = 0.509. At v'/v' c < 2.0, 
i.e. in the homogeneous phase, a peak is found at 9 w 0. This 
is consistent with densely grafted polymers in the homoge- 
neous phase in good solvents, in which the grafted polymers 
stretch away from the substrate due to the excluded volume in- 
teraction between adjacent polymers [13, 28]. In contrast, the 
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regions of 2-dimensional phase separation, i.e. v'/v' c > 3.0, 
reveal another sharp peak at 6/n * 0.4, in which the peak at 
=s is covered. This indicates that the grafted polymers tend 
to be stretched laterally, as expected. This result corroborates 
the location of the binodal line given in Fig. 2. 
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FIG. 5. /(£>) atp' p /p' pc = 0.509 and 1 x 10 6 MCS for various values 
of v'/v' c . f{0) = for 6/n > 0.5, due to the hard substrate placed on 
the xy-plane. 

In conclusion, single-component linear homopolymers ho- 
mogeneously grafted onto a planar substrate have been sim- 
ulated in 3 dimensions using the molecular Monte Carlo 
method. The 2-dimensional microphase separation of grafted 
polymers in the lateral direction of the substrate has been ex- 
amined and a variety of lateral domain patterns have been re- 
vealed. This result is similar to that for the 3-dimensional mi- 
crophase separation of non-grafted diblock copolymer melts. 
The structural phase diagram of a homopolymer brush has 
been constructed, the binodal line of which is qualitatively 
consistent with that for the non-grafted homopolymer system. 

Both non-grafted block copolymer melts and grafted ho- 
mopolymers exhibit microphase separation and various do- 
main patterns. A block copolymer is composed of chemi- 
cally immiscible blocks that are separated into microstruc- 
tures. However, such phenomena are also observed for the 
grafted homopolymers, despite the absence of chemical im- 
miscibility along each polymer chain and distinct molecular 
architectures. This indicates that the topological constraints, 
i.e. the grafting points of homopolymer brushes and the con- 
necting points between the blocks in each block copolymer, 
are the origin of microphase separation and the various do- 
main patterns. This also demonstrates that both intramolecu- 
lar topological constraints, i.e. the connecting points between 
the blocks, and the extramolecular topological constraints, i.e. 
the grafting points, yield these phenomena. This result is a 
universal characteristic of arbitrary polymers. 
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